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Abstract. We consider the problem of computing all-pairs shortest paths 
in a directed graph with real weights assigned to vertices. 
For snnxiiO-1 matrix C, let Kc be the complete weighted graph 
on the rows of C where the weight of an edge between two rows is equal 
to their Hamming distance. Let MWT(C) be the weight of a minimum 
weight spanning tree of Kc- 

We show that the all-pairs shortest path problem for a directed graph 
G on n vertices with nonnegative real weights and adjacency matrix Ac 
can be solved by a combinatorial randomized algorithm in tim^fl 



0(nJn + mm{MWT(A G ), MFFT(A' G )}) 

As a corollary, we conclude that the transitive closure of a directed graph 
G can be computed by a combinatorial randomized algorithm in the 
aforementioned time. 

We also conclude that the all-pairs shortest path problem for uniform 
disk graphs, with nonnegative real vertex weights, induced by point sets 
of bounded density within a unit square can be solved in time 0(n 2 ' 75 ). 



1 Introduction 

The problems of finding shortest paths and determining their lengths are funda- 
mental in algorithms. They have been extensively studied in algorithmic graph 
theory. A central open question in this area is if there is a substantially subpu- 
bic in the number of vertices algorithm for the all-pairs shortest path problem 
for directed graphs with real edge weights (APSP) in the addition-comparison 
model [24 27 . For several special cases of weights and/or graphs substantially 
subcubic algorithms for the APSP problem are known [3 8 23 25 26 27]. How- 
ever, in the general case the fastest known algorithm due to Chan [8] (see also 



3 The notation 0( ) suppresses polylogarithmic factors and B l stands for the trans- 
posed matrix B. 



[5]) runs in time 0(n 3 log 3 log n/ log 2 n), achieving solely a moderate polylog- 
arithmic improvement over the 0(n 3 ) bound yielded by Floyd- Warshall and 
Johnson's algorithms |1I27] . 

The situation is different for directed graphs with real vertex weights. Re- 
cently, Chan has shown that the APSP problem for the aforementioned graphs 
can be solved in time (3(n 2 - 844 ) [5] and Yuster has slightly improved the latter 
bound to 0(n 2 - 842 ) by using an improved bound on rectangular multiplication 

The basic tool in achieving substantially subcubic upper bounds on the run- 
ning time for the APSP for directed graphs with constrained edge weights or 
real vertex weights are the fast algorithms for arithmetic square and rectangu- 
lar matrix multiplication |10I14] . One typically exploits here the close relation- 
ship between the APSP problem and the so called distance or (min, +) product 
13123124125127125] . 

Unfortunately, these fast algorithms for matrix multiplication, yielding equally 
fast algorithms for Boolean matrix product, are based on recursive algebraic ap- 
proaches over a ring difficult to implement. Thus, another central question in 
this area is whether or not there is a substantially subcubic combinatorial (i.e., 
not relaying on ring algebra) algorithm for the Boolean product of two n x n 
Boolean matrices |4|22 '"24i . Again, the fastest known combinatorial algorithm 
for Boolean matrix product due to Bansal and Williams [1] running in time 
0(n 3 log 2 log n/ log 9 '' 4 n) achieves solely a moderate polylogarithmic improve- 
ment over the trivial 0(n 3 ) bound. On the other hand, several special cases of 
Boolean matrix product admit substantially subcubic combinatorial algorithms 
|5ll3l2()j . 

In particular, Bjdrklund et al. 5, provided a combinatorial randomized al- 
gorithm for Boolean matrix product which is substantially subcubic in case the 
rows of the first n x n matrix or the columns of the second one are highly clus- 
tered, i.e., their minimum spanning tree in the Hamming metric has low cost. 
More exactly, their algorithm runs in time 0(n(n + c)), where c is the minimum 
of the costs of the minimum spanning trees for the rows and the columns, re- 
spectively, in the Hamming metric. It relies on the fast Monte Carlo methods 
for computing an approximate minimum spanning tree in the L\ and L2 metrics 
given in |16|17j . 

The assumption that the input directed graph is highly clustered in the 
sense that the minimum spanning tree of the rows or columns of its adjacency 
matrix in the Hamming metric has a subquadratic cost does not yield any direct 
applications of the algorithm of Bjorklund et al. [5] to shortest path problems, 
not even to the transitive closure. The reason is that the cost of the analogous 
minimum spanning tree can grow dramatically in the power graphs of the 
input graph. In particular, we cannot obtain directly an upper time-bound on 
the transitive closure of Boolean matrix corresponding to that for the Boolean 
matrix product from [5] by applying the asymptotic equality between the time 



4 In the i-th power graph there is an edge from v to u if there is a path composed of 
at most i edges from v to u in the input graph. 



complexity of matrix product over a closed semi-ring and that of its transitive 
closure over the semi- ring due to Munro [3T]. The reason is the dependence of 
the upper bound from [5] on the cost of the minimum spanning tree. 

In this paper, we extend the idea of the method from [5] to include a mixed 
product of a real matrix with a Boolean one. We combine the aforementioned 
extension with the ideas used in the design of subcubic algorithms for important 
variants of the APSP problem 3 26 , in particular those for directed graphs with 
vertex weights |8l25j . to obtain not only a substantially subcubic combinatorial 
algorithm for the transitive closure but also for the APSP problem in highly 
clustered directed graphs with real vertex weights. 

For an n x n — 1 matrix C, let Kq be the complete weighted graph 
on the rows of C where the weight of an edge between two rows is equal to 
their Hamming distance. Let MWT{C) be the weight of a minimum weight 
spanning tree of Kq- We show that the all-pairs shortest path problem for 
a directed graph G on n vertices with nonnegative real weights and an ad- 
jacency matrix Aq can be solved by a combinatorial randomized algorithm in 
0(n 2 ^n + mm{MWT(A G ), MWT(A* G )}) time. It follows in particular that the 
transitive closure of a directed graph G can be computed by a combinatorial ran- 
domized algorithm in the aforementioned time. 

Our algorithms are of Monte Carlo type and by increasing the polylogarith- 
mic factor at the time bounds, the probability that they return a correct output 
within the bounds can be amplified to 1 ^r, where a > 1. 

Since there are no practical or combinatorial substantially subcubic-time algo- 
rithms not only for the APSP problem but even for the transitive closure problem 
for arbitrary directed graphs at present, our simple adaptive method might be a 
potentially interesting alternative for a number of graph classes. 

As an example of an application of our method, we consider the APSP prob- 
lem for uniform disk graphs, with nonnegative real vertex weights, induced by 
point sets of bounded density within a unit square. We obtain a combinatorial 
algorithm for this problem running in time 0(y^rn 2 ' 75 ), where r is the radius of 
the disks around the vertices in a unit square. 

The recent interest in disk graphs, in particular uniform disk graphs, stems 
from their applications in wireless networks. In this context, the restriction to 
point sets of bounded density is quite natural. In Fiirer and Kasiviswanathan 
provided a roughly 0(rt 2 5 )-time preprocessing for approximate 0(y / n)-time dis- 
tance queries in arbitrary disk graphs. 

Our paper is structured as follows. In the next section, we show a reduction 
of the APSP problem for directed graphs with real vertex-weights to a mixed 
matrix product of a distance matrix over reals with the — 1 adjacency matrix. 
In Section 3, we present an algorithm for such a mixed product which generalizes 
that for the Boolean matrix product from [5] and runs in subcubic time if the 
input — 1 matrix is highly clustered. By combining the results of Sections 2,3, 
we can derive our main results in Section 4. In the next section, we present 
the application of our method to uniform disk graphs induced by point sets of 
bounded density. We conclude with final remarks. 



2 A Reduction of APSP to Mixed Matrix Products 



2.1 The APSP problem 

Formally, the All-Pairs Shortest Paths problem (APSP) in a directed graph 
G = (V, E) with real weights w(v) associated to vertices v € V is to compute 
the |V| x |V| distance matrix D G such that Dg{v,u) is the distance Sg{v,u) 
from v to u in G, i.e., the minimum total weight of vertices on a path from v to 
u in G. An additional goal of the APSP problem is to compute a concise data 
structure representing the shortest paths. 

Note that Sg(v, it) is equal to the minimum total weight of inner vertices on 
a path from v to u in G increased by the weights of v and u. 

We shall assume \V\ — n throughout the paper. 

For i = 0,1,..., n — 1, let 6 G (v,u) be the distance from v to u on paths 
consisting of at most i edges, i.e., the minimum total weight of vertices on a 
path from v to u having at most i edges in G. Next, let D l G be the |V| x \V\ 
matrix such that D l G [v,u] is equal to 5 G (v,u). 

For convention, we assume S G (v, v) = and S G (v, u) = +oo for v ^ u. Hence, 
D G has zeros on the diagonal and +oo otherwise. In D G , all the entries D G [v, u] 
where (v, u) E E are set to w(v) + w(u) instead of +oo. Thus, both D G and D G 
can be easily computed in time 0(n 2 ). 

2.2 Mixed Matrix Products 

Let A be an n x n matrix over R U {+oo}, and let B be an n x n matrix with 
entries in {0, 1}. The mixed right product C of A and B is defined by 

C[i,j] = mm{A[i,k]\l <k<nk B[k,j] = 1} U {+00} 

If C[i, j] +00 then the index k such that C[i, j) = A[i, k] (and thus B[k, j] — 
1) is called a witness for C[i,j}. Analogously, the mixed left product C' of B and 
A is defined by 

C'[i,j] = mm{A[k,j]\l <k<nk B[i,k] = 1} U {+00}, 

and if C'[i,j] 7^ +00 then the index k such that C'[i,j) = A[k,j] is called a 
witness for C'[i,j}. 

An n x n matrix W such that whenever C[i,j] 7^ +00 then W[i,j] is a 
witnesses for C[i,j] is called a witness matrix for the right mixed product of A 
and B. Analogously, we define a witness matrix for the left mixed product of B 
and A. 

2.3 The Reduction 

Let A G denote the n x n adjacency matrix of G = (V,E), i.e., A G [v,u] = 1 iff 
(v,u) e E. 



Lemma 1. For an arbitrary i e {0, 1, — 2}, D l G x can be computed on the 
basis of D l G and the right mixed product of D G with Aq or D G and the left mixed 
product of Aq with D G in time 0(n 2 ). 

Proof. It is sufficient to observe that for any pair v, u of vertices in G, Dq' 1 [v, u] 
is equal to 

mm{D G [v, u],mm{D G [v,x] + w(u)\l <x<nk A G [x,u] = 1} U {+oo}} 
Symmetrically, D G ~ 1 [v,u] is equal to 

mm{D G [v,u],min{D G [x,u) + w(v)\l < x < n k A G [v, x] = 1} U {+oo}} 

□ 

The following lemma follows the general strategy used to prove Theorem 3.4 
in@. 

Lemma 2. Let G be a directed graph G on n vertices with nonnegative real 
vertex weights. Suppose that the right (or left) mixed product of an n x n matrix 
over RU {+00} with the adjacency matrix Aq of G along with the witness matrix 
can be computed in time T m i X (n) — Q(n 2 ). The APSP problem for G can be 
solved in time (^{n 1 ^ yjT m i x {n)). 

Proof. We begin by computing D^ 1 for some t G [2, n] which will be specified 
later. By Lemma [1] this computation takes time 0(tT m i X (n)). 

It remains to determine distances between pairs of vertices where any shortest 
path consists of at least t edges. For this purpose, we determine a subset B of 
V, the so called bridging set [26], hitting all the aforementioned long paths. We 
apply the following fact to I = t and sets of t vertices on shortest consisting of 
exactly t—1 edges, similarly as in [3=8-2 5126] . 

Fact 1. Given a collection of N subsets o/{l, where each subset has size 

exactly I, we can find a subset B of size 0{{n/l) log n) that hits all subsets in the 
collection in time O(Nl). 

Since our application of Fact 1 is analogous to those in [318125126] . we solely 
sketch it referring the reader for details to the aforementioned papers. 

Note that for each pair v, u, of vertices for which any shortest path has at 
least t edges there is a pair v',u' of vertices on a shortest path from v to u 
such that any shortest path from v' to u' has exactly t — 1 edges. For all such 
pairs v',vf, we can find a shortest path on t — 1 edges, and thus on t vertices, 
by backtracking on the computation of -D^T 1 and using witnesses for the mixed 
products. In total, we generate 0{n 2 ) such paths on t vertices in time 0{tn 2 ). 
The application of Fact 1 also takes time 0(tn 2 ). 

Next, we run Dijkstra's single-source shortest path algorithm pQ for all ver- 
tices in the bridging set B in the input graph G and in the graph resulting from 
reversing the direction of edges in G. In this way, we determine Dg[v,u] for all 
pairs (v, u) e (B x V) U (V x B). 



Now, it is sufficient for all remaining pairs (v, u) in V x V to set 
Dq[v, u] = imn{D t G ~ 1 (v, u), min{D<3[?;, b] + Dc[b, it] — w(b)} 

in order to determine the whole F>g- 

The computation of D 1 ^ 1 takes 0(tT m i X (n)) time which asymptotically is 
not less than the 0(tn 2 ) time taken by the construction of the bridging set. The 
runs of Dijkstra's algorithm and the final computation of Dg require 0(jn 2 ) 
time. By setting t — T " 3 ( T { j ; we obtain the lemma. □ 

3 Fast Computation of the Mixed Products for Clustered 
Data 

Our algorithm for the right (or, left) mixed product relies on computation of an 
approximate minimum spanning tree of the columns (or rows, respectively) of 
the Boolean input matrix in the Hamming metric. 

3.1 Approximate Minimum Spanning Tree in High Dimensional 
Space 

For c > 1 and a finite set S of points in a metric space, a c- approximate minimum 
spanning tree for S is a spanning tree in the complete weighted graph on S, with 
edge weights equal to the distances between the endpoints, whose total weight 
is at most c times the minimum. 

In [16] (section 4.3) and [15] (section 3), Indyk and Motwani in particular 
considered the bichromatic e-approximate closest pair problem for n points in R 
with integer coordinates in O(l) under the L p metric, p £ {1,2}. They showed 
that there is a dynamic data structure for this problem which supports insertions, 
deletions and queries in time 0(dn 1 ^ 1+e ^) and requires 0(dn + n 1+1 /( 1+c ))- 
time preprocessing. In consequence, by a simulation of Kruskal's algorithm they 
deduced the following fact. 

Fact 2. For e > 0, a 1 + e-approximate minimum spanning tree for a set of n 
points in R d with integer coordinates in 0(1) under the L\ or metric can be 
computed by a Monte Carlo algorithm in time 0(dn 1+1 ^ 1+e ^). 

In |17j Indyk, Schmidt and Thorup reported even slightly more efficient (by a 
poly- log factor) reduction of the problem of finding a 1+ e-approximate minimum 
spanning tree to the bichromatic e-approximate closest pair problem via an easy 
simulation of Prim's algorithm. 

Note that the L\ metric for points in R n with 0, 1-coordinates coincides 
with the n-dimensional Hamming metric. Hence, Fact 2 immediately yields the 
following corollary. 

Corollary 1. For e > 0, a 1 + e-approximate minimum spanning tree for a set 
o/nO-1 strings of length n under the Hamming metric can be computed by a 
Monte Carlo algorithm in time 0{n 2+1 ^ 1+ ^). 



3.2 The Algorithm for Mixed Matrix Product 



The idea of our combinatorial algorithm for the right mixed product C of A 
with B and its witness matrix is a generalization of that from [5]. Let P(r,v) 
denote a priority queue (implemented as a heap) on the entries A[r, k] such that 
B[k,v] = 1 ordered by their values in nondecreasing order. 

First, we compute an approximate minimum spanning tree of the columns of 
B in the Hamming metric. Then, we fix a traversal of the tree. Next, for each row 
r of A, we traverse the tree, construct P(r, start) where start is the first column 
of B in the tree traversal and then maintain P(r, v) for the currently traversed v 
by updating P(r, u) where u is the predecessor of v in the traversal. A minimum 
element in P(r, v) yields a witness for C[r, v}. The cost of the updates in a single 
traversal of the tree is proportional to the cost of the tree modulo a logarithmic 
factor. 



f Algorithm 1 \ 

Input: n x n matrix A over flu{+oo} and an n x n Boolean matrix B; 
Output: A witness matrix W for the right mixed product C of A and B. 
Comment: P(r,v) stands for a priority queue on the entries A[r, k] s.t. B[k, v] = 1 
ordered by their values in nondecreasing order. 

1. Compute an 0(logn)-approximate minimum spanning tree Tb of the columns of 
B in the Hamming metric; 

2. Fix a traversal of the tree Tb linear in its size; 

3. Set start to the first node of the traversal; 

4. For each pair of consecutive neighboring columns v, u in the traversal, precompute 
the set D v<u of positions where Is occur in v but not in u and the set D u<v of 
positions where Is occur in u but not in v; 

5. For each row r of A do 

— Construct the priority queue P(r, start) and if P(r, start) 7^ set W[r, start] 
to the index k where A[r, k] is the minimum element in P(r, start); 

— Traverse the tree Tb and for each node v different from start compute the 
priority queue P(r,v) from the priority queue P(r,u), where u is the prede- 
cessor of v in the traversal, by utilizing D VtU and D UiV . If P(r,v) 7^ set 

\^ W[r, v] to the index k where A[r, k] is the minimum element in P(r,v). ^ / 

Lemma 3. Algorithm 1 is correct, i.e., it outputs the witnesses matrix for the 
right mixed product of matrices A and B. 

For an n x n Boolean matrix C, let Kc be the complete weighted graph on 
the rows of C where the weight of an edge between two rows is equal to their 
Hamming distance. Next, let MWT(C) be the weight of a minimum weight 
spanning tree of Kc- 

Lemma 4. Algorithm 1 can be implemented in time 0(n(n+MWT(B t )))+t(n), 
where t(n) is the time taken by the construction of the O (log n)- approximate 
minimum weight spanning tree in step 1. 



Proof. Step 1 can be implemented in time t(n) while steps 2,3 take time 0(n). 
Step 4 takes 0(n 2 ) time. The block in Step 5 is iterated n times. 

The first step in the block, i.e., the construction of P(r, start) takes 0(n log n) 
time. The update of P(r,u) to P(r,v) takes 0(\ogn(\D ViU \ + \D UiV \)) time. Note 
that |-Dt,. u | + |-Du,-u| is precisely the Hamming distance between the columns v 
and u. It follows by the O(logn) approximation factor of Tg that the total time 
taken by these updates is O^MWTiB 1 ) log 2 n). 

We conclude that Step 5 can be implemented in time 0(nMWT(B t )). □ 

Theorem 1. The right mixed product of two n x n matrices A over RU {+00} 
and B over {0, 1} can be computed by a combinatorial randomized algorithm in 
time 0{n{n+MWT{B t ))). Analogously, the left mixed product of B and A can be 
computed by a combinatorial randomized algorithm in time 0(n(n+MWT(B))). 

Proof. By Corollary [TJ an 6>(log n)-approximate minimum spanning tree can be 
constructed by a Monte Carlo algorithm in time 0(n 2 ) (observe that n 1 ^ = 0(1) 
if / = J?(logn)). Hence, by Lemmata 02 21 we obtain the theorem for the right 
mixed product. The upper bound on the time required to compute the left mixed 
product follows symmetrically. □ 

4 Main results 

Lemma [5] combined with Theorem [T] yield our main result. 

Theorem 2. Let G a directed graph G onn vertices with nonnegative real vertex 
weights. The all-pairs shortest path problem for G can be solved by a combinato- 
rial randomized algorithm in time 0(n 2 ^Jn + min{MWT (Aq) , MWT(A t G )}). 

By setting vertex weights, say, to zero, we obtain immediately the following 
corollary. 

Corollary 2. The transitive closure of a directed graph G on n vertices can be 
computed by a combinatorial randomized algorithm in time 
0(n 2 + mm{MWT(A G ), MWT(A t G )}) . 

Equivalently, we can formulate Corollary [5] as follows. 

Corollary 3. The transitive closure of an n x n Boolean matrix B (over the 
Boolean semi-ring) can be computed by a combinatorial randomized algorithm in 
time 0{n 2 yjn + min{MWT(B), MWT(B t )}). 

5 APSP in vertex- weighted uniform disk graphs of 
bounded density 

In this section, we consider uniform disk graphs that are induced by a set P of 
n points in a unit square in the plane that are 6(n)-dense, where b : N —> N. 



Formally, we say that P is 6(n)-dense iff each cell of the regular sjn x sjn grid 
within the unit square contains at most b(n) points. The vertices of such an 
induced disk graph are the points in P, and two vertices are adjacent in the 
graph iff their Euclidean distance is at most r, where r is a positive constant not 
exceeding 1. We shall term the aforementioned graphs as uniform disk graphs 
induced by b(n) -dense point sets. 

Lemma 5. Given two intersecting disks on the plane of the same radius r with 
the distance d between centers, the area of the symmetric difference is 0(rd). 




Proof. AC = ^,AB = r. 

The area of the triangle ABC is 

Area ABC = l -ACBC = ~yj - J = ±dy/4r*-<P. 

The area of the circular sector ABD is 

I i 

Area A BD = -r 2 ZBAC = -r 2 arccos 

The area of BCD is AreascD = AreaABD — AreaABC- 
The area of the symmetric difference 

Area = 2(nr 2 — AAreascD) = 27rr 2 — 4r 2 arccos (j^J + 4r 2 — d 2 . 

Finally, by using Taylor series expansion 

**»»»(s)-**(f-5 + ((s)"))- 

= 2nr 2 - 2dr + 0(d 2 ) = 2?rr 2 - 2dr + 0{rd) 
and V 4r 2 — d 2 < 2r we get Area = 0(rd). □ 




Lemma 6. Let G be a uniform disk graph induced by a b(n)-dense point set. 
For each edge (v,u) of G, the number of vertices in G that are a neighbor of 
exactly one of the vertices v, u, i.e., the Hamming distance between the two 
rows in the adjacency matrix of G corresponding to v and u, respectively, is 
0(r x b(n) (dist(v, u) x n + y/n)). 

Proof. The number of vertices of G that are a neighbor of exactly one of the 
vertices v and u is at most the minimum number of cells of the regular y/n x y/n 
grid within the unit square that cover the symmetric difference S(v, u) between 
the disks centered at v and u, respectively, multiplied by bin). The aforemen- 
tioned number of cells is easily seen to be at most the area A(v,u) of S(v,u) 
divided by the area of the grid cell, i.e., A(v,u) x n, plus the number of cells of 
the grid intersected by the perimeter of S(v, u), i.e., 0(ry/n). By Lemma[3J we 
have A(v, u) — 0(dist(v,u) x r). Hence, the aforementioned number of cells is 
0(r(dist(v, u) x n + y/n)). 

□ 

The following lemma is a folklore (e.g., it follows directly from the upper 
bound on the length of closed path through a set of points in a d-dimensional 
cube given in Lemma 2 in [18 ). 

Lemma 7. The minimum Euclidean spanning tree of any set of n points in a 
unit square in the plane has total length 0{\/n). 

Combining Lemmata [51 El we obtain the following one. 

Lemma 8. For a uniform disk graph G induced by a b(n) -dense n-point set, 
a spanning tree of the rows ( or, columns ) of the adjacency matrix of G in the 
Hamming metric having cost 0(rn 3 ^ 2 ) can be found in time 0(n 2 ). 

Proof. Construct a minimum Euclidean spanning tree of the n points forming 
the vertex set of G. It takes time O(nlogn) and the resulting tree T has total 
length 0(y/n) by Lemma [7] Form a spanning tree U of the rows (or, columns) 
of the adjacency matrix of G by connecting by edge the rows corresponding to 
v and u iff (v, u) E T. By LemmalHland the 0(y/n) length of T, the total cost of 
U is 0{rn 3 / 2 b(n)). □ 

By plugging Lemma [8] into Theorem [21 we obtain our main result in this 
section. 

Theorem 3. Let G be a uniform disk graph, with nonnegative real vertex weights, 
induced by a b(n) -dense n-point set. The all-pairs shortest path problem for G 
can be solved by a combinatorial algorithm in time 0(\/rn 2 y/b(n)). 

In the application of the method of Theorem [5] yielding Theorem [31 we can 
use the deterministic algorithm of Lemma[5]to find a spanning tree of the rows or 
columns of the adjacency matrix of G instead of the randomized approximation 
algorithm from Fact 2. 



By straightforward calculations, our upper time-bound for APSP in vertex- 
weighted uniform disk graphs induced by 0(l)-dense point sets subsumes that for 
APSP in sparse graphs based on Dijkstra's single-source shortest-path algorithm, 
running in time 0(nm), where m is the number of edges, for r » nT 1 ^ . 

Finally, we can also easily extend Theorem [3] to include uniform ball graphs 
in a d-dimensional Euclidean space. In the extension, the term y/r in the upper 
time-bound generalizes to v H -1 . 

6 Final Remarks 

We can easily extend our main result to include solving the APSP problem for 
vertex and edge weighted directed graphs in which the number of different edge 
weights is bounded, say by q. This can be simply achieved by decomposing the 
adjacency matrix Aq into the union of up to q matrices A\, A2, ...Ai in one-to- 
one correspondence with the distinct edge weights and consequently replacing 
each mixed product with I such products in Lemmata [I] [2] In the final up- 
per bound, MWT{A G ) and MHT(4) are replaced by £)j =1 MWT(Ai) and 

Ei=i MWT{A\), respectively 

It is an interesting problem to determine if there are other natural graph 
classes where MWT(Ag) or MWT(A f G ) are substantially subquadratic in the 
number of vertices. 

It follows from the existence of the so called Hadamard matrices [7J that 
there is an infinite sequence of graphs with rii x m adjacency matrices A4 such 
that mm{MWT(Ai),MWT(A\)} = n{{ ni ) 2 ) holds. 
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